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ABSTRACT 

We study the simulated CO emission from elliptical galaxies formed in the mergers of gas-rich disk 
galaxies. The cold gas not consumed in the merger-driven starburst quickly resettles into a disk-like 
configuration. By analyzing a variety of arbitrary merger orbits that produce a range of fast to slow- 
rotating remnants, we find that molecular disk formation is a fairly common consequence of gas-rich 
galaxy mergers. Hence, if a molecular disk is observed in an early- type merger remnant, it is likely the 
result of a "wet merger" rather than a "dry merger". We compare the physical properties from our 
simulated disks (e.g. size and mass) and find reasonably good agreement with recent observations. 
Finally, we discuss the detect ability of these disks as an aid to future observations. 
Subject headings: Galaxies: Elliptical and Lenticular, cD — Galaxies: evolution — Galaxies: Kinematics 
and Dynamics — ISM: molecules 
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1. INTRODUCTION 

The "red and dead" nature of early-type galaxies 
such as ellipticals and lenticulars has long been part of 
the standard lore. The general lack of cold molecular 
gas conducive to the formation of stars in these galax- 
ies (|Lees et al.lll991h is typically considered the reason 
for their redder colours (paldry et al.ll2QQ3 ). Recently, 
however, a number of studies have shown that some 
molecular gas can exist in early- types, o f tentimes in a 
disk-like configuration (ITnoue et al.l [19961: iRupenl 1 19971: 
Wiklind et ainmliWound 12002b lOkuda et all [20051 : 



Das et aLll2Q0a : lYoundl2005l : lYoung et al.ll2008f ). 



A wide variety of mechanisms have been proposed 
for the origin of molecular disks in early- type galaixes. 
These include external origins s uch as gas accreted 
from the intergalactic medium (Khochfar & B urker 
1 20051) or gas remnant after a merger (|Lees et al.l 1 19911 
IKnapp fc Rupenlll996l ). Some theories have additionally 
speculated an int ernal origin, i.e. from the mass loss 
of ev olved stars (Fa ber fc Gallagherl 119761 : iCiotti et al.l 
Il991f ). Since star formation occurs in the presence 
of molecular gas, the properties of these disks can be 
used to study the weak star formation that occurs in 
as many as 30% of local early-type galaxies (|Yi et al.l 
120051 : iKavirai etHI I2007D . Their kinematics and dy- 
namics can also be used to model and infer the po- 
tentials of their ho st galaxies (|de Zeeuw fc Franxl 119891 : 
iCretton et alll2000f ). 

Given that ellipticals are thought to form in the merger 
of two gal axies where grav itational disruption is preva- 
lent (e.g. ICox et aTll2006L and references therein), the 
idea that either stellar or molecular disks may exist in 
elliptical galaxies seems counterintuitive. There is, how- 
ever, some theoretical evidence to suggest that stellar 
disks are born after galaxy mergers (iBarnes fc Hernauistl 
1996; Roberts onet al.l 120061 : iRoberston fc Bullockl 120081 



Hopkins et al. 2009). Since stars form from molecular 
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gas, it is conceivable that molecular disks may also form 
after "wet" (gas-rich) galaxy mergers (as opposed to 
"dry" or gas-poor mergers). 

In this paper we use smoothed particle hydrodynamic 
(SPH) simulations of gas-rich major mergers combined 
with molecular CO line radiative transfer calculations in 
order to investigate the potential formation of molecu- 
lar disks in early- type galaxies. This approach allows 
us to compare synthetic CO observations of our model 
galaxies, which span the full range between fast and 
slow-rotators, with data from the literature. Although 
major mergers are rare events in the context of LCDM 
(iFakhouri fc Mall2008l : iGuo fc : Whitdl2fifi3 IStewartetall 
120081 ). the simulations presented here can still be viewed 
as useful limiting cases. Our paper is organized as fol- 
lows: In $21 we describe our numerical methods; In $3] 
and we discuss the formation of molecular disks in 
our simulations and compare these with recent observa- 
tional results, and in $5] we conclude. 

2. SIMULATIONS AND MODELS 

2.1. Simulations 

We model major mergers between two disk galaxies 
using a modified version of the publicly available N- 
body/SPH code gadget-2 (|Springel2QQ5h . Detailed de- 
scriptions of the parent galaxies as well as the multiphase 
IS M, star formation and blackhole fee dbac k can be found 
injSpringel fc Hernquistl (|2002L l2QQ3h and lSpringel et al.l 
(|2005af ). We refer the reader to these works for a more 
detailed description of the simulations. The most rel- 
evant aspects to this work will be summarized briefly 
below. 

The parent disks in our simulations h ave dark mat- 
ter halos initialized with iHernquistl (jl990f ) profiles. They 
have a concentration index of c = 9, a spin parameter 
of A = 0.033 and a circular velocity of ^200 = 160km/s. 
A total of 120,000 dark matter particles and 80,000 disk 
particles were used to model each progenitor galaxy. The 
disk consists of 40% gas by mass, with the remainder be- 
ing collisionless star particles. The total (halo and bary- 
onic) mass of a parent galaxy is ~ 2 x 1O 12 M . Softening 
lengths of 100 pc and 200 pc were used for baryons and 
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dark matter respectively. 

The key evolutionary steps during the major merger 
simulations are as follows. After the initial passage of the 
merging gal axies, tidal torquing drives gas in to the cen- 
tral regions ([Barnes fc Hernqu ist 1991. fl996[ ), leading to 
a period of enhanced star formation rates. As the galax- 
ies undergo final coalescence during the second passage, 
the merger-driven sta rburst drives the gal axy through a 
LIRG/ULIRG phase ([Younger et al.ll2QQ9f ). A combina- 
tion of gas consumption and AGN feedback quench the 
star formation, and the final merg er product evolves pas - 
sively into a red elliptical galaxy ([Springel et al.ll2005b[ ). 
We analyze the results of the simulations for four identi- 
cal mergers which vary only in orbit. The merger orbit 
angles are chosen arbitrarily, and are given in Table [U 

2.2. Radiative Transfer 

In order to investigate the simulated CO properties 
of the remnants of our model mergers, we post-process 
the SPH output using the 3D non-L TE molecular line 
radiative tr ansfer code, TURTLEBEACH ([Narayanan et all 

HooS HooH). 

The molecular gas fraction in galaxies is determined 
by various different properties such as metallicity, dust 
content, interstellar radiation fields, density and tem - 
perature ([Hollenbach et al.lll971l : [Pelupessv et al.ll2QQ6[ ). 
However, since the spatial resolution is limited in our 
simulations, this cannot be calculated explicitly as the lo- 
cations of individual stars and gas clouds are unknown. 
Hence, we assume that half of the cold neutral gas in 
each grid cell is atomic and the other half molecular. 
This is consis tent with st u dies o f nearby star-forming 
galaxies (e.g. iKeres et al.l ([20 03)). We also examined 
a sub-sample of our simulations by scaling the H2 gas 
fraction by the ambient interstellar pressure accor ding 
to local scaling relations (Blitz & Ro solowskvl l2QQ6h . In 
the simulations examined here, the differences in the two 
methods for assigning the molecular gas mass are negli- 
gible. 

The molecular abundance leve ls are set uniformly at 
Galactic values ([Lee et al.l Il996[ ). The relatively high 
abundance of C O in the galaxy (CO/H 2 = 1.5 x 10 -4 , 
iLee et al.l (JT996)) in comparison to other H2 tracers such 
as HCN, CS and HCO + , makes it the best observational 
tracer for molecular hydrogen. Hence, we focus on the 
radiative transfer of CO, particularly the readily observ- 
able CO (J=l-0) transition in this work. 

TURTLEBEACH calculates the molecular level popula- 
tions, and hence the source functions by considering 
both radiative excitations/deexcitations as well as col- 
lisions with H2. In practice, a solution grid to the level 
populations is guessed and model photons are emitted 
isotropically via a Monte Carlo sampling. By assuming 
statistical equilibrium and balancing radiative and colli- 
sional excitation, de-excitation and stimulated emission, 
the level populations are updated until convergence is 
reached. The source function can be subsequently de- 
termined. The emergent intensity is then calculated by 
integrating the equation of radiative transfer through the 
model grid. For f urther details on the un derlying algo- 
rithms, please see Na rayanan et al.l (|2008[ ). 

3. RESULTS 




Fig. 1. — Synthetic JHK colours mapped onto RGB for two 
arbitrary viewing angles of our fiducial simulated galaxy, model 
d3e. The images were produced with the publicly available dust 
radiative transfer code, Sunrise. The variable input parameters 
included a stellar birthcloud covering fraction of /pdr = 0.3 and 
an evolving dust mass with the metallicity following a dust-to- 
metals ratio of 0.4. Evidently, the model galaxy has an elliptical 
morphology. 



To determine when the product of the major merger 
has evolved into an elliptical, we examine the star forma- 
tion rates (SFRs) and stellar morphologies of the simu- 
lated system. We pick snapshots well after the merger- 
induced starburst, when the stars have relaxed into a 
spheroid-type configuration to serve as our model early- 
type galaxies. The SFRs are typically of order <0.5 M 
yr _1 . To confirm the early- type morphology of our sim- 
ulated galaxies, we have run dust radiative transfer cal- 
culations in post-processing to determine the synthetic 
SEDfl We plot the simulated JHK colours mapped 
onto RGB for model d3e (hereafter our fiducial mo de0) 
at two arbitrary viewing angles in Figure [TJ A study by 
ICox et all J2006) using the same simulations as this work 
shows that the early-type galaxies formed span the full 
range between fast and slow-rotators. Their rotation pa- 
rameters (Vmaj/cr, where V ma j is the major axis rotation 
speed and Vmag/cr > 1 defines a fast-rotator) range be- 
tween 0.06 to 1.03 when derived from a model including 
dissipation. 

In Figure O we show centroid velocity maps in CO 
(J=l-0) for one sight-line (where the molecular disk is 
highly inclined) through our fiducial model and a sim- 
ilarly chosen sight-line through model d3f. The maps 
are 8 kpc on a side. The velocity-integrated intensity is 
shown for comparison. 

We search for evidence of molecular disks in the final 
elliptical by examining the CO centroid velocity maps 
from 100 random viewing angles. Ordered rotation, as 
evidenced by clear blue and red centroid velocity peaks 
when a disk is highly inclined, is observed along several of 
the sight- lines for each merger orbit angle implying that 
the molecular gas is in fact located within a rotating disk. 
The fact that this is observed for all of our merger orbit 
angles indicates that molecular disk formation may be a 
relatively common occurence in merger-induced elliptical 
galaxy formation. 

The possibility of molecular disks being a common oc- 
currence in elliptical galaxies is interesting. First, early- 

4 For this, we employ the publicly avai lable code SUNRISE; a 
full description of the code can be found in Jons son et al.l {2010). 
For the variable input parameters, we chose a stellar birthcloud 
covering fraction of /pdr = 0.3 and a dust-to-metals ratio of 0.4. 

5 For reference with previous pa pers in the litera t ure, t his is 
additionally the fidu cial model in Naravanan et al. ( 2008]) and 
lYounger et a l. (2009). 
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Fig. 2. — (top) Centroid velocity maps in CO (J=l-0) for a nearly 
edge-on sight-line through the final elliptical galaxy in our fiducial 
model (d3e) and d3f. The units on the velocity are km/s. (bottom) 
Same as top, except these are velocity-integrated intensity maps 
instead. The units on the velocity-integrated intensity are km/s 
• erg/s/cm 2 /Hz/steradian. The clear distinction between the red 
and blue velocity peaks in the centroid velocity plot indicates the 
presence of a disk. As we find evidence for molecular disk formation 
in the ellipticals produced through all of our merger orbit angles, 
we conclude that a variety of merger configurations are conducive 
to disk formation. 



type galaxies are canonicaly thought to be relatively de- 
void of molecular gas. Second, if early-types form in hi- 
erarchical merging of dynamically cold disks, one might 
naively expect the disks to be disrupted in the merging 
process. 

The molecular disks in our simulations form from gas 
not consumed in the merger-driven starburst. During 
the merger, gravitational torques on the gas (by the 
stars) cause the gas to lose anuglar momentum and 
collapse toward the cen ter ([Barnes fe Hernquistl 119961 : 
iMihos fc Hernau ist 1996). A large increase in central gas 
density ensues leading to a burst of star formation. The 
remnant gas, originally in the outer regions of the pro- 
genitors, re-virializes into a rotationally-supported disk 
around the compact central potential induced by the 
post-merger stars. A similar formation mechanism for 
post-mer ger stellar disk fo r matio n has been discussed in 
detail by Robertson et al. (120061) . iRoberston fc Bullockl 
(|2008D and lHopkins et al.l (|2009l ). 

4. DISCUSSION 
4.1. Comparison with Observations 

lYoung et all (j2008|) (hereafter Y08) observed the CO 
(J=l-0) emission in four early- type galaxies using the 
BIMA array, which has a primary beam size of 100" . Of 
the four, three showed evidence for molecular disks. The 
properties of these three galaxies are listed in Table [2j 
We compare our results to their observation by placing 
our model galaxies at 17 Mpc (~ the mean distance of 
the Y08 galaxies) and convolving the CO (J=l-0) model 
spectra with a 100" circular Gaussian beam. We list the 
relevant properties derived from our models in the same 
table for comparison. 

We define the effective radius, R ei to be the radius 
at which the stellar mass profile drops from its cen- 



tral value by a factor of 1/e. This implicitly assumes 
a constant mass- to- lig ht ratio at these smaller radii 
([Cappellari et al.ll2006f ). From comparing values in Ta- 
ble^ we find the effective radius of NGC 3032 is similar 
to that of the d3e merger-remnant. The effective radius 
of NGC 4459 is in between that of the d3e and d3k rem- 
nants. The effective radius of NGC 4526 is similar to 
that of the d3k remnant. 

The size of the molecular disks in our simulations is 
determined in the following way. For a near edge-on 
sight-line for each merger orbit angle, a line is put down 
connecting the maximum and minimum velocities for the 
molecular disks in the final ellipticals. An interpolation is 
performed to find the centroid velocities along this line 
at various distances from the center of the disk (taken 
to be where the velocity equals zero along the line). The 
resulting rotation curve of the molecular disk is sinusoidal 
in shape as would be expected for a rotating disk. We 
take the radius corresponding to the extremum of the 
rotation curve as the radius of the disk rdisk- These 
values are listed in Table [T] and alongside the Y08 values 
in Table [2j Reasonable correspondence is seen between 
the observed values in Y08 and our models to within a 
factor of ~2-3. 

Finally, the simulated molecular gas masses within the 
disks are compared (see Tabled]) to those observed. The 
study of Y08 uses a CO-to-H2 conversion factor of 3x 10 20 
cm -2 (K km/s) -1 when deriving their molecular gas 
masses. In comparing these molecular gas masses (see 
Table [2]), the agreement seems poorer at first glance. The 
typical discrepancies are within a factor of ~4-10. How- 
ever, this is about equal to the dispersion seen in the CO- 
to-H2 conversion factor in local galaxies (which ranges by 
a factor of 2-5 within the Gala xy and by up to an order of 
magnitude in other galaxies; [ Downes & Solomon 1998; 
Gao fc Soloinonl 120041 : 1 Arimoto et aH ll996: Boselli et al.l 
2002ft . 

The total (baryonic + dark matter) mass of the 
merger-remnant for each model is 3.9 x 1O 12 M . The 
percentage of the total mass enclosed within the molec- 
ular disk radius is given under the column heading 
in Tableland are typically of order ~l-2%. This result 
may be applicable in determining the total mass of an 
observed galaxy. 

4.2. Detectability 

Finally we investigate the detectability of these molec- 
ular disks using two different approaches: from their syn- 
thetic CO spectra and from their centroid velocity maps. 

From the 100 arbitrary model viewing angles, we esti- 
mate the detectability of the disk as the fraction of these 
sight-lines where an obvious broadened and double-side- 
peaked line (due to rotation) is visible. An automated 
peak finding routine is used to test whether each spec- 
trum has this characteristic multi-peak profile. If it does, 
then we consider the disk to be detectable. 

Examples of such spectra obtained by placing the 
galaxy at varying distances and convolving with different 
beam sizes are shown in Figure [3] with 20 km/s smooth- 
ing (a typical resolution used by Y08). The green line in 
the plot corresponds to the models we compared to the 
observations of Y08 in §4.11 (i.e. distance to galaxy =17 
Mpc, beam size = 100"). The blue and red lines repre- 
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Fig. 3. — CO (J=l-0) spectra for a sight-line through the molec- 
ular disk of model d3m. The blue line is the expected spectrum of 
a galaxy 17 Mpc away convolved with a 50" beam. The green line 
is the expected spectrum of a galaxy 17 Mpc away convolved with 
a 100" beam. The red line is the expected spectrum of a galaxy 34 
Mpc away convolved with a 100" beam. The characteristic signa- 
ture of a disk is apparent in the broadened lines with two smaller 
peaks on opposite sides indicating rotation. We label these disks as 
"detectable" and define the detectability of the disk as the fraction 
of sight-lines through the disk that have these suggestive profiles. 
Since the CO line profile may change due to beam size or distance 
to the galaxy, we vary these two parameters in the plot to test the 
robustness of this characteristic line profile. 



sent galaxy distances and beam sizes of 17 Mpc, 50" and 
34 Mpc, 100" respectively. Although this is an arbitrary 
sight-line through the molecular disk of model d3m, it is 
representative of all the models in this study. 

The detectability of disks using spectral line profiles 
for the models comparable to Y08 are given in Table [1] 
under the heading 'DETdisk (spectra)'. The sightline av- 
eraged detectability of disks (e.g. the fraction of viewing 
angles over which the characterstic 3-peaked line profile 
is visible) is ^60-75%. This is comparable to the de- 
tection fraction of molecular disks in early-types by Y08 
(75%). The smaller beam size and more distant galaxy 
cases are detectable to similar percentage levels (within 
-2-5%). 

The spectroscopic CO line profile may change in shape 
due to beam size and distance to a galaxy. For example, 
if the galaxy is isolated and unresolved in a telescope 
beam, as the beam size gets larger the line tempera- 
ture will decrease while the line shape stays the same. 
If the galaxy is resolved, then as the beam gets bigger 
or smaller, more or less velocity components within the 
galaxy are included or excluded so the line shape may 
change. Figure [3] and the similar detectabilities suggests 
that this technique for detecting underlying molecular 
disks is relatively robust against a variety of modeled 
galaxy distances and beam sizes. 

Another method that can be used to detect molecu- 
lar disks is by analyzing the centroid velocity data. In 
this case, we can estimate the detectability of the disk 
as the fraction of sight-lines in each merger model that 
show a single distinct blue velocity peak and a single 



distinct red velocity peak, like in Figure [2j If the cen- 
troid velocity map has these distinct peaks, we consider 
the disk detectable. Another automated peak finder is 
used to perform this search. We find similar detectabil- 
ities to the spectra method for the models comparable 
to Y08. These are given in Table [1] under the heading 
L DET disk (CVy and mostly lie in the —60-75% range as 
before. 

5. CONCLUSIONS 

We study the early- type galaxies formed in gadget- 2 
simulations of gas-rich major mergers at four arbitrary 
merger orbit angles. A 3D non-LTE radiative transfer 
code creates 100 CO (J=l-0) spectra for each merger 
remnant by "viewing" it at arbitrary viewing angles. 
Molecular gas disks are detected in all four merger rem- 
nants which range between fast to slow-rotators. This 
implies that molecular gas disks, thought to reform 
after a gas-rich merger through a re-virialization pro- 
cess induced by the compact central potential of post- 
merger/post-starburst stars, may not be a rare occur- 
rence in early- type galaxies. Consequently, if an early- 
type merger remnant contains a molecular disk, it is 
likely that a wet merger was responsible for its forma- 
tion rather than a dry merger. 

The rare nature of major mergers in the standard 
LCDM context though, means that the simulations stud- 
ied here are likely not perfect analogs to most early-type 
galaxies in the local universe. They are, however, inter- 
esting limiting cases in the formation of molecular disks 
in early- type galaxies. We would expect that a more 
typical merger ratio like 1:3 will leave a n SO- like rem- 
nant with < 20% of the mass in the bulge ([Hopkins et al.l 
l2009f ). This implies a smaller disruption to the original 
gas disks and hence, a higher probability of having a 
molecular disk in the remnant. 

We compare the molecular disk radii and the molec- 
ular gas masses enclo sed within t h ese ra dii to the ob- 
servational results of I Young et al.l (|2008f ). Reasonably 
good agreement between our models and their results 
are found to within a factor of ~2-3 for disk radius. The 
molecular gas masses differ by a factor of ~4-10. How- 
ever, these differences are within the uncertainty of the 
CO-to-H2 conversion factor. 

We find the total mass enclosed within the molecular 
disk radius to be ~l-2% of the total mass (baryons + 
dark matter) of the galaxy. This result may be useful for 
estimating the total mass of an observed galaxy. 

We investigate the detectability of underlying molecu- 
lar disks using two methods: from CO spectra and cen- 
troid velocity maps. From model CO spectra, we derive 
typical disk detection rates of ~60-75%. This compares 
well to the 75% det ection rate of molec ular disks in the 
galaxies observed by I Young et al.l ([20081 ) . Using the cen- 
troid velocity map approach, a similar disk detection rate 
is achieved. 



The computations in this paper were run on the 
Odyssey cluster supported by the FAS Sciences Division 
Research Computing Group at Harvard University. 
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TABLE 1 

Simulation orbit angles and derived properties of product ellipticals 



Model 


Orbit Angle 

4>1 01 4>2 #2 


T(iisk 

(kpc) 


M(r < r disk ) 
(M@) 


M H2 (r < r disk ) 
(M ) 


/ 

(%) 


DETdisk 
(spectra) 


DETdisk 
(CV) 





d3e 


30 


60 


-30 


45 


2.0 


4.1 x 10 10 


5.3 x 10 7 


1.1 


0.47 


0.50 


d3m 








71 


90 


2.0 


3.7 x 10 10 


2.1 x 10 9 


0.95 


0.61 


0.66 


d3f 


60 


60 


150 





3.6 


7.2 x 10 10 


3.2 x 10 9 


1.9 


0.72 


0.66 


d3k 


-109 


30 


71 


30 


1.8 


3.7 x 10 10 


2.3 x 10 9 


0.95 


0.71 


0.75 



Note. — Models are referred to by the designations under the "Model" heading throughout the 
paper. 



TABLE 2 

Comparison of properties between Y08 galaxies and our model galaxies 



Property 


NGC 3032 


NGC 4459 


NGC 4526 


d3e 


d3m 


d3f 


d3k 




Distance (Mpc) 
R e (kpc) 
M H2 (r < r d isk) (M Q ) 
rdisk (kpc) 


21.4 
0.93 
5.0 x 10 8 
1.5 


15.7 
2.7 
1.6 x 10 8 
0.67 


16.4 
3.5 
5.7 x 10 8 
1.1 


17 
1.2 
5.3 x 10 7 
2.0 


17 
6.0 
2.1 x 10 9 
2.0 


17 

4.8 
3.2 x 10 9 
3.6 


17 
3.6 
2.3 x 10 9 
1.8 



